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RESIDUAL HUBBLE-BUBBLE EFFECTS ON SUPERNOVA COSMOLOGY 
Benjamin Sinclair 1 , Tamara M. Davis 1 ' 2 , Troels Haugb0lle 3 

ABSTRACT 

Even in a universe that is homogeneous on large scales, local density fluctuations can imprint a 
systematic signature on the cosmological inferences we make from distant sources. One example 
is the effect of a local under-density on supernova cosmology. Also known as a Hubble-bubble, it 
has been suggested that a large enough under-density could account for the supernova magnitudc- 
rcdshift relation without the need for dark energy or acceleration. Although the size and depth of 
under-density required for such an extreme result is extremely unlikely to be a random fluctuation 
in an on-average homogeneous universe, even a small under-density can leave residual effects on our 
cosmological inferences. It is these small under-densities we consider here. 

In this paper we show that there remain systematic shifts in our cosmological parameter measure- 
ments, even after excluding local supernovae that are likely to be within any small Hubble-bubble. 
We study theoretically the low-redshift cutoff typically imposed by supernova cosmology analyses, 
and show that a low-redshift cut of zq ~ 0.02 may be too low based on the observed inhomogeneity 
in our local universe. 

Neglecting to impose any low-redshift cutoff can have a significant effect on the cosmological pa- 
rameters derived from supernova data. A slight local under-density, just 30% under-dense with scale 
70/i _1 Mpc, causes an error in the inferred cosmological constant density J7a of ~ 4%. Imposing a 
low-redshift cutoff reduces this systematic error but does not remove it entirely. A residual systematic 
shift of 0.99% remains in the inferred value £7a even when neglecting all data within the currently pre- 
ferred low-redshift cutoff of 0.02. Given current measurement uncertainties this shift is not negligible, 
and will need to be accounted for when future measurements yield higher precision. 
Subject headings: cosmology: observations, theory — supernovae: general 



1. INTRODUCTION 

To estimate global properties of the universe, such as 
the average matter density, we make models that include 
those properties as parameters and test which parame- 
ter values best fit the data. The simplest cosmological 
model that provides a good fit to supernova data is the 
ACDM model, in which the universe is assumed to be 
homogeneous and dominated by cold dark matter and a 
cosmological constant, parametrized by their normalized 
densities f^M and The assumption of homogene- 

ity however, is only valid over large scales. Evidently, 
on the scales of stars, galaxies, and clusters of galaxies 
the matter distribution is far from homogeneous. Such 
inhomogeneities give rise to varying expansion rates in 
different parts of the universe, and varying gravitational 
potentials. 

Thus, as light makes its journey from source to ob- 
server it experiences not a cosmological redshift due to 
a single rate of expansion H{t), but rather varying de- 
grees of cosmological redshift determined by the local 
rate of expansion H (r, t). In addition, the spatially vary- 
ing gravitational potential contributes differing gravita- 
tional red- and blue-shifts, as the light travels through 
under- and over-dense regions. 

To a large extent these effects cancel each other out 
over the path of the photon. In a flat universe containing 
only matter the canceling is exact, however in the pres- 
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ence of other contributions, such as from dark energy or 
curvature, the potential wells evolve as the light passes 
through them in such a way so as to alter the total red- 
shift. This effect is known as the Integrated Sachs- Wolfe 
(ISW) effect and has been recently used to detect the 
presence of dark energy from temperature enhancements 
in the CMB corr elated with the density of galaxies along 
the line of sight ([Giannant onio et al.H2008h . 

Nevertheless the ISW effect is small (an order of mag- 
nitude below the level of the primordial fluctuations in 
the CMB) and as we look over a range of directions we 
benefit from spatial averaging which further smooths any 
variation due to density fluctuations along any particular 
line of sight. Spatial averaging also compensates for the 
over- or under-density at the sources. 

The one irreducible effect comes from any density fluc- 
tuation inside which we happen to reside. Then the ab- 
solute depth of the potential (rather than just the change 
in depth during a light-crossing time) comes into full ef- 
fect. In that case the photon has no chance to climb back 
out of the potential well and feels the full brunt of the 
density fluctuation. Whether we be in an under-dense 
or over-dense region, it is clear that the fluctuation light 
travels through immediately before reaching us is not 
compensated for, as it would have been were the light 
left to continue its journey. 

This final redshift will add a systematic shift in the 
magnitudc-rcdshift data to which we fit our homogeneous 
models, and may impede our ability to deduce the aver- 
age density of the universe outside our local bubble. 

We are therefore motivated to study the possible effect 
of a local inhomogeneity on our cosmological inferences. 
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In this paper we focus on supernova cosmology. We gen- 
erate simulated data with a local under-density then fit it 
with an homogeneous cosmological model to test whether 
the derived parameters match the input parameters for 
the average density outside the local void (Sect.HJ). Su- 
pernova cosmologists typically reject low-redshift data on 
the basis that peculiar velocities are relatively high com- 
pared to the cosmological recession velocities and Hubble 
bubble effects most prominent. We therefore implement 
a low-redshift cutoff, denoted Zq, and investigate how 
well the fit converges to the input model as this low- 
redshift cutoff changes (Sect. [5]). 

Initially we use a large void in a matter-only universe 
to demonstrate the concept clearly, then in Sect. [472] in- 
troduce a cosmological constant as well. In Sect. 14.31 we 
assess the likely impact of voids of the size we expect to 
find in random fluctuations in a standard ACDM cosmol- 
ogy. Our conclusions are analyzed in Sect. [71 We now 
begin by first discussing the nature and sizes of voids in 
our universe. 

2. VOID SIZES 

The word 'void' has been used in the literature to mean 
several different things. The 'voids' we refer to here are 
gaussian under-densities, of the kind you might expect 
to naturally occur as random fluctuations in a gener- 
ally homogeneous universe. While very large voids are 
non-gaussian in nature, as a first approximation, and for 
small voids, the assumption of gaussianity should hold. 
Our motivation is to study the possible effect of a local 
density fluctuation of a size that is likely to occur in a 
standard ACDM cosmology. We choose to concentrate 
on under-densities rather than over-densities because of 
the interest they have already generated in the commu- 
nity: firstly due to the possibility of mimicking a cosmo- 
logical constant (for very large scale under-densities) and 
secondly because of the observational indications from 
supernova cosmology and galaxy surveys that there may 
be a local void (a 'Hubble bubble'). 

The ev idence for a local Hubble bubble was first 
found by IZehavi et al.l (|1998|) . in which the supernova 
data appeared to show that Hq within 70/i -1 Mpc was 
6.5%±2.2% higher than the value outside that distance 
(presuming a flat 17m =1.0 model). They explained this 
by a 20% under-density surrounded by a dense shell, and 
noted that this size corresponds rou ghly to the size ob - 
served for local large scal e structure (IGeller et al!l ll997l. 
such as the 'great walls'. Uha. Riess. fe Kirshnerl (2007) 
refined this observation with more supernovae, and con- 
sidered a void embedde d in a ACDM univer se. They find 
the significance of the IZehavi et al.l (|1998l) result drops 
to S H = (H in - H out )/H out ~4.5%±2.1% in the ACDM 
case and find a similar value of 6.5%±1.8% with their 
new data. 

These resu lts have since been c hallenged by several pa- 
pers, such as lConlev et al.l (|2007l ) who showed that using 
a different light-curve-fitter with a different color treat- 
ment can remove the evidence of a void. Data from 
other distance-measurement techniques do not find a 
sig nificant Hubble b ubble . Tully-Fisher measurements 
by iGiovanelli et al.l ()1999l ) find a statisti cally insignifi- 
cant H ubble bubble of Sh ~1%±2% while iHudson et ah! 
(2004) use the peculiar velocity field within 120/i _1 Mpc 
to find modest evidence of locally enhanced Hubble ex- 



pansion with Sh ~ 3 ± 1.3%. 

These observations have all been measurements of the 
local expansion rate. Alternatively one can look to large 
galaxy redshift surveys and measure directly from the 
observed stru cture the typical s i ze of v oids in the galaxy 
distribution. iHovle fe Vogelev l (12004? ) found using the 
2dFGRS that 35% of the universe consists of voids larger 
than 10/i _1 Mpc, with a mean galaxy under-density of 
<5 g ai = (pin - Pout)/pout = -0.94 ± 0.02. Fluctuations in 
the dark matter density will be less extreme than this, 
since this is a measurement of the galaxy density fluctu- 
ation, which is enhanced due to galaxy bias (galaxies are 
at the denser parts of the dark matter distribution) . 

Perhaps the most relevant observations are those of our 
local universe. In a su rvey of our local un iverse, out to a 
redshift of z ~ 0.15, IGeller et al.1 (|1997t ) find evidence 
for inhomogeneity on the 100/i _1 Mpc scale. Earlier 
IGeller fe Huchral (|1989tl noted that voids of 50/i _1 Mpc 
with density only ~ 20% of the mean are ubiquitous 
to all surveys, and the largest local structure, the 'Great 
Wall' has dimensions at least 60/i _1 Mpc by 170ft -1 Mpc. 

Theoretical calculation s of the void size that is likely 
in a ACDM universe (e.g. iFurlanetto fe Pir"anll2006D find 
void sizes and density contrasts somewhat smaller than 
these and than the 2dFGRS survey observations. This 
remains a point of tension. Moreover, recent observa- 
tions of bulk flows in the universe, such as those by 
iKashlinskv et all ()2008l ) using the Sunyaev-Zeldovich ef- 
fect in the CMB, show that large-scale motions also ap- 
pear to be larger than predicted by theoretical ACDM 
calculations and simulations. 

This tension between observation and theory, married 
with the hint of a Hubble bubble in the supernova data 
and the observed distance to the 'great walls' of our lo- 
cal structure, suggests that a void as large as 70/i _1 Mpc 
is a reasonable, and interesting, void size to test. The 
wide range of density contrasts considered in the litera- 
ture makes it difficult to define a 'typical' under-density. 
The depth of under-density depends strongly on the size 
of void considered. We choose a value of 6 — —0.3 (cor- 
responding to pi n = 0.7/Oout) for all the simulated voids 
in this paper. The typical peculiar velocity correspond- 
ing to such a void is of the order of 100km s _1 . Taking 
this as a monopole velocity at a scale of 70/i _1 Mpc it 
is a typical (a ~1-1. 5 sigma level) fluctuati on in a stan- 
dard ACDM model (jHaugbolle et all 120071) . and having 
the additional coincidence of being near the Centre of the 
void with correlated monopole velocities out to at least 
70 h~ 1 Mpc does not lower this probability much more 
than to the 2 sigma level. 

The precise value of the size and density contrast of 
the void we test is not of prime importance. We focus 
on the qualitative results which remain roughly the same 
regardless of choice of void size. In particular we concen- 
trate on how a void of any size affects supernova cosmol- 
ogy and whether we can remove its effect by neglecting 
data within the void. 

An alternativ e cosmologica l model has recently 
emerged (see e.g. lEnqvistl 120081 and references therein), 
where we are at the center at a gigaparsec sized void 
embedded in an Einstein-de Sitter universe. It has 
been shown that while we have to give up the Coper- 
nican Principle, and fine tuning is needed to have 
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Figure 1. A snapshot at to of the radial profile of the matter den- 
sity f2ivi(f, to) used to generate the simulated data sets in Sects. |4~TI 
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Void models used. 
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the observer very near the center of the void, such 
a model is viable when compared to current state- 
of-the-art cosmological observation s (lAlnes et al.1 120061: 
' Garcia-Bellido fe Haugballel l2008atfbL 120091: iZibin et al.1 
20081 iClifton et all 12003 lAlexander et all 120091 ). That 
includes not only supernova, but also CMB an d BAO 
observations (|Garci'a-Bellido fc Haugb0llg I2009D . and a 
void of such dimensions, would be difficult to spot in 
galaxy redshift surveys, due to the small gradient in 
the density profile. The model is attractive because 
there is no dark energy and the observed late-time ac- 
celeration is a consequence of a larger Hubble rate near 
the center of the local under- density. The main draw- 
back is the required coincidence of having the Milky 
Way near the center to make the CMB radiation from 
the surface of last scatterin g appear close to isotropic 
(EI nes fc Am arzguiou MM)- 

We begin by considering a large void of scale size 
ro = 700ft.~ 1 Mpc to demonstrate the concept clearly. A 
void of this order of magnitude is what would be required 
to explai n the supernova data without a cosmological 
constant (| Alexander et all 12009). We embed this large 
void in both a flat f2 m = 1.0 universe (without a cosmo- 
logical constant, Sect. 14. and in a ACDM universe with 
(Cl M ,Q A ) = (0.3,0.7) (Sect. S3). We assess the likely 
impact of a more realistic sized void of ro = 70/i~ 1 Mpc 
in Sect. 14.31 A summary of the void sizes we consider 
appears in Table [TJ and is depicted in figure [TJ 



3. INVESTIGATING THE HOMOGENEOUS UNIVERSE IN 
THE PRESENCE OF A LOCAL UNDER-DENSITY 

In this section we explain the method we use to mea- 
sure the effect of a local inhomogeneity on the cosmo- 
logical parameters derived using supernova cosmology. 
We generate simulated data containing an under-dense 



region at the origin, and then investigate the impact of 
neglecting to account for the under-density by fitting the 
data with an homogeneous ACDM model. The impact 
of the void is the difference between the best-fit cosmo- 
logical parameters derived using the homogeneous model 
and those input into the LTB simulation. 

3.1. Modeling a local under-density 

The model we use to describe a local void is a Lemaitrc- 
Tolman-Bondi (LTB) model, which is isotropic but inho- 
mogeneous. The metric of our model can be written 

(a'r + a) 2 



ds z 



-dt z 



1 — k(r)r' 



-dr z + a 2 r 2 dQ 



(1) 



where a(r, t) is an effective scale factor, k(r) describes 
the curvature as a function of coordinate distance r, and 
prime denotes differentiation with respect to r. 

Defining the transverse and longitudinal Hubble pa- 
rameters H{r, i) = f and H L (r, t) = ^r±| from the Ein- 
stein equatio ns one can construct an effective Fri edmann 
equation (see IGarcia- Bellido fc Haugb0lld[2OO8ah 



H 2 (r,t)=H 2 (r) 



n M (r) 



ao(r) 



a(r,t) 

+ (1 - fi M (r) - fi A (r)) 



(2) 



a(r,t) 



+ ^A(r) 



where Hq{t) = H(r,to), and ao(r) = a(r,to) = 1 is 
a gauge freedom giving the scale factor at to- The to- 
tal matter density Qj^ (r) and the cosmological constant 
£Ia(t) are related to Ho(r), the curvature k(r), and the 
physical matter density pm at t = to as, 



8ttG 



H 2 (r)aUr)r* 
8itG 

" A(r) = 3i^ A ' 
k(r) =H 2 (r)(n M (r) 



dr'r' 2 p M (r\t ), (3) 



•n A (r)-l)a§(r). 



(4) 
(5) 



The time to big bang as a function of distance is a func- 
tion of radius, which can be seen by integrating ([2]) with 
respect to a(r,t), but only models with a constant time 
to big bang give a well motivated gr owing void profile as 
a function of time (|Zibin et al. 2008) , and below we only 
consider such models. 

Using the scale factor we can determine the time of 
emission for light at a given redshift by solving the red- 
shift equation, 



dt 



1 



d(log(l + z)) H L {r{z),t{z)Y 



(6) 



We can also determine the luminosity distance (and 
hence apparent magnitude) from 



d L (z) = (1 + z) 2 a[r(z),t(z)}r(z). 



(7) 



The distance modulus is given by p = 51og 10 (dL(.z)) + 25, 
with luminosity distance measured in Mpc. 

In our model the matter distribution at large distances 
is a constant f^*, but nearby the density decreases to 
fl 1 ^ at the origin. Thus, the LTB model reproduces the 



4 



Sinclair, Davis, Haugb0lle 



important features we are investigating: an on-average 
homogeneous universe in which we reside within a local 
Hubble-bubble. The LTB model can accept any density 
profile for the local inhomogeneity. We choose to imple- 
ment a simple gaussian, according to, 

n M (r)=ns? t + - n^)e- r2 / r o, (8) 

= n^(l+de- r2 / r o), (9) 

where 6 = (fift - figf*)/^ = (p™ - is the 

density contrast of the void. Given the density profile 
the Hubble parameter now Ho(r) is uniquely defined 
up to a proportionality, d ue to the requirement of a 
consta nt time to big-bang ([Garcia- Bellido & Haugbeilld 
l2008a| K We consider models in which the homogeneous 
universe outside the void consists exclusively of matter 
(fijj* = 1.0) as well as models in which the universe ex- 
terior to the void is homogeneous ACDM , (fi^ 4 , fi A ut ) = 
(0.3,0.7). 

In the latter case the mass density distribution is still 
given by equation and the asymptotic cosmological 
constant density is 

n A (r -> oo) = 0.7. (10) 

Close to the center of the void, the Hubble expansion rate 
is higher, while the value of the cosmological constant, 
Pa, is by construction constant, giving lower values for 

n A (r). 

3.2. Fitting an homogeneous model in the presence of 
an under- density 

The redshift and distance modulus calculated from 
equations [6] and [7] are the output from the LTB model, 
which we treat as the data we would collect if we were 
indeed situated at the center of a Hubble-bubble in an 
otherwise homogeneous (over large scales) universe. Ex- 
plicitly, we input a set of z values, which we take to be the 
redshifts of a set of supernovae we are observing, numer- 
ically integrate along a null geodesic in the LTB space- 
time to obtain the position and time of emission from the 
source at redshift z, and then use this position and time 
to calculate a distance modulus, from equations [2] and 
To solve the equati ons we use an extended version of th e 
easyLTB program ([Garcia- Bellido fc Haugbollc 2008a), 
that can handle an arbitrary mixture of dark matter and 
cosmological constant in the void. 

For the sake of finding the best fit cosmological pa- 
rameters, we do not add simulated observational errors 
to the 'correct' distance modulus, the reason for which 
is explained at the end of this section. To perform the 
\ 2 minimizing fit, we simply assign an uncertainty to 
the distance modulus of ±0.2 mag, which is a typical 
uncertainty in the measurement of distance modulus for 
supernovae. 

We then fit a set of parametrized homogeneous models 
to the LTB data and find which ACDM cosmology is the 
best fit (by minimizing \ 2 )- This process is entirely anal- 
ogous to the methodology cosmologists use to measure 
cosmological parameters from real data. The homoge- 
neous models are parametrized by the values of (OM,OL, 
w) and we marginalize over HO as per the standa rd pro- 
cedure in supernova cosmology (Kim ct al. 20041). 

If the local inhomogeneity had no effect on our deduc- 
tions then we would find that the best fit set of Q.m and 



Qa would match the fi^vT* and f2 A ut we had input into 
our LTB model. 

When we come to consider goodness-of-fit of the best 
fit models in section 16.11 it is necessary to generate a 
set of observational errors to obtain realistic values of 
goodness-of-fit. We generate observational errors accord- 
ing to a normal distribution with standard deviation 0.2 
mag and add them to the distance modulus. We re- 
peat the procedure for different sets of normally dis- 
tributed errors and average the result to get an mean 
goodness of fit, mean best fit matter density, mean best 
fit cosmological constant, and mean best fit equation of 
state. Note, that the errors are normally distributed, 
and therefore the mean best fit values of each parameter 
should converge (with infinite sets of different normally 
distributed errors) to the value obtained when no nor- 
mally distributed errors were added. This is the reason 
why we only need to add the observational errors when 
we want to estimate the goodness-of-fit; and why in sec- 
tions 14.11 to 14.31 we can avoid running multiple fits by 
performing a single fit without simulated errors, but giv- 
ing each data point a 0.2 mag observational uncertainty. 

3.3. Sensitivity of Result to Size of Data Set 

It is important to note that the measure we use for 
goodness-of-fit will go down rapidly for the homogeneous 
model as the number of data points is increased and/or 
the precision of each data point is improved. This is 
because an incorrect model (the homogeneous model in 
this case) can not mimic the correct model in the limit 
of infinite data (or vanishing uncertainty j3- 

This also means that the best fit parameters change 
as the number of data points varies. When a model is 
a good description of the data, the data points will be 
evenly distributed above and below the theoretical curve. 
However, when the model does not reflect the data there 
will be regions in which all data points lie above (or be- 
low) the theoretical curve. (See the low-redshift regions 
of Fig. [2]for an example.) Increasing the number of data 
points therefore gives more weight to the regions in which 
the model is the worst fit, since these regions have the 
most significant impact on x 2 ■ The best fit theoretical 
curve becomes distorted towards the data in the poor- 
fitting regions, at the expense of moving further away 
from the data in regions where the theoretical curve was 
a better fit. 

Thus, when the test model is not an accurate repre- 
sentation of the system, as in our case of fitting a ACDM 
model to LTB data, the cosmological parameters we in- 
fer will change depending on how many data points we 
choose to use. 

We have chosen to use 301 data points spread out uni- 
formly with redshift over the redshift range < z < 1.7, 
since that is approximately t he number ava i lable to cur- 
rent supernova surveys (e.g. iKessler et al.l [200 9) Fl As 
more supernova data are collected, and our uncertainty 
in their measurements reduces, it will become increas- 

4 As long as the model does not include the correct model in 
some limit of its parameters. 

5 We could also choose to weight the redshift distribution ac- 
cording to the number of supernovae currently observed at each 
redshift, but that is unnecessary given the general nature of the 
analysis here. 
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Figure 2. In this Hubble diagram distance modulus relative to the 
empty universe is plotted against redshift. This shows the ACDM 
model (dashed line) that is the best fit to the magnitude-redshift 
data generated from the LTB model (dotted line) with ro = 0.7 
h~ 1 Gpc and Hi*! = 0.7 as used in Sect. I4.ll The data was given 
an uncertainty of 0.2mag. 

ingly apparent whether an inhomogeneous or homoge- 
neous model gives the better fit. 

4. RESULTS: NEGLECTING TO ACCOUNT FOR A LOCAL 
UNDER-DENSITY 

4.1. Large void in matter-only universe 

We begin by setting the void size to an extremely large 
ro = 700/i _1 Mpc. This is far larger than the expected 
scale of inhomogeneities in the standard ACDM model, 
but is approximately the minimum size required to ex- 
plain the supernova data without invoking a cosmological 
constant. We use it here to clearly demonstrate the ef- 
fect an under-dense region can have on our cosmological 
inferences. 

We fit the standard ACDM cosmology to our simulated 
void data and find the best fit to be (^m, ^a) = (0.83 ± 
0.08,0.39). This model is plotted against the LTB data 
in Figure [5] 

The best fit matter density, = 0.83, is significantly 
lower than the actual value of f^ 1 * = 1.0 for the exter- 
nal universe in the LTB model. This simplified example 
demonstrates that a Hubble-bubble creates the illusion 
of a lower matter density than the true density outside 
the void. 

The best fit model also has a significant cosmological 
constant, J7a = 0.39, whereas the input model had none, 
Q A ut = 0.0. In section IfTTl we demonstrate that this ex- 
tra parameter improves the goodness of fit significantly 
compared to a model with J7a = 0.0. This demonstrates 
that the addition of an extra parameter (the cosmolog- 
ical constant) will be strongly supported when falsely 
attempting to fit an homogeneous model to data with 
such a large local void0 

4.2. Large void in ACDM universe 

Having demonstrated the systematic shift that occurs 
if observations are made from a region of Dm lower than 

6 Note that we chose a void of this size specifically because of this 
feature. To recover the ACDM parameters (^Mi^a) = (0.3,0.7) 
some fine tuning is required on the shape of the void. The shape 
must be made sharper than our gaussian profile allows. 
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Figure 3. This hubble diagram shows the ACDM model (dashed 
line) that is the best fit to the magnitude-redshift data generated 
from the LTB model containing a cosmological constant (dotted 
line) with r = 0.7 /i^Gpc and (njj.fij™*) = (0.21,0.3) as used 
in Sect.rOl 



the surroundings, we now show that the same behavior 
is expected in a universe with a cosmological constant. 

The model we now consider consists of an under-dense 
region embedded in the currently preferred ACDM uni- 
verse with (Om, ^a) = (0.3, 0.7), which is a situation we 
may indeed find ourselves in if ACDM is the correct de- 
scription of the universe. However, this model is still 
not 'realistic' in the sense that we still use a large void 
size, tq =700/i _1 Mpc, which is much larger than is likely 
to occur in a ACDM universe, so the model is not self- 
consistent. We use this model for comparison with the 
large void considered in the previous section. 

Figure [3] depicts the best fit ACDM model, which has 
(r>M,^A) = (0.28,0.86). This corresponds to a system- 
atic error of 6.7% and 23% respectively. This example 
differs from the pure matter model tested in the previous 
section, because the homogeneous model we are fitting 
to the data has the same number of parameters as the 
simulated data exterior to the local void. We find that 
in this case a non-zero A is even more important to the 
fit than when there was no external A in the simulation. 
In Table we show the results for the best fits of homo- 
geneous modelfl with and without A. In each case we 
show the x 2 values, indicating the goodness of fit. In all 
cases the \ 2 1S much worse when no cosmological con- 
stant is included (as expected since adding parameters 
always allows a better fit), but the improvement in fit 
is more dramatic when the simulated universe included 
A. That is as expected since in this case A actually does 
relate to something real in the model. In Sect. 16.11 we 
study in more detail how much the additional parame- 
ters improve the fits, and whether the improvements are 
enough to fool us into believing false parameters. 

The obvious question arises as to whether a local void 
can dupe us into believing in a dark energy equation of 
state that differs from the cosmological constant value of 
w = — 1. When we consider the more realistic sized voids 
in the next section we also fit a dark energy model in 
which w is allowed to vary, wCDM. This again restores 

7 Note, in Table [2] we have used a different fitting methodol- 
ogy where we simulate observational errors, giving slightly differ- 
ent best fit parameters. The difference is within the estimated 
uncertainties of our best fit parameters. 
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Small A 
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0.25 
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301-3 


1.001 




Small A 


Flat tcCDM 


0.30 


0.70 


-1.1 


298.8 


301-2 


0.999 





us to the position of having one extra parameter that 
may allow a better fit. In essence we are exchanging 
the ro parameter, which describes the size of the void in 
our input LTB model, with a false parameter, a free w. 
The extra parameter is guaranteed to allow a better fit 
to the data. The crucial point for us is whether the fit 
is improved so much that the extra parameter appears 
justified. If so, then a local void is a significant danger 
for misleading our cosmological inferences. 



4.3. Small void in ACDM 



universe 



Having demonstrated that the Hubble-bubble alters 
the inferred external matter density and cosmological 
constant density, we now investigate the magnitude of 
this effect for a more modest sized under-density, of the 
kind that is predicted in self-consistent models of struc- 
ture formation in ACDM and seen in the typical size of 
structures actually observed in our universe. The exact 
size of typical under-densities in our universe is the topic 
of ongoing debate, which we summarized in Sect. [U We 
have chosen to test a void with r =70/i _1 Mpc, and den- 
sity contrast 5 = —0.3, as representative of the 'likely' 
under-density we may find ourselves in. This corresponds 
to having a maximal monopole velocity of 120 km s _1 in- 
side the void, which is a typical monopole velocity for a 
shell in a standard ACDM model (|Haugb0lle et ali l2007). 

Imprinting this void on a flat , fl° A ut ) = (0.3,0.7) 
background, we find that the best fit ACDM model has 
parameters (Om,^a)=(0.299, 0.73). The error in the 
best fit parameters is vastly reduced for the realistic sized 
under-density model compared to the large void model 
examined in Sect. 14.21 but is still large enough to be of 
concern. The errors in the best fit matter density and 
best fit cosmological constant density are 0.3% and 4.3% 
respectively. These errors must be evaluated in the con- 
text of the current uncertainty in the cosmological pa- 
rameters. Using a combination of CMB, BAO and SN 
data lKomatsu et al.l ((2001I report n M h 2 = 0.1358±£'oo36 
and fl A = 0.726 ± 0.015 (Table 1, WMAP+BAO+SN 
mean values). These correspond to errors of 2.7% and 
2% respectively. 

Compared to the observational error in the cosmolog- 
ical constant density, the error incurred by fitting an in- 
adequate model can be significantly larger. The corre- 
sponding systematic error in the matter density is small 
but not negligible compared to the current uncertainty, 
and as measurements become more precise the relative 
importance of this potential systematic error will in- 
crease. In Sect. [5] we discuss how this error is mitigated, 
but not removed, by introducing a low-z cutoff on the 
data. 
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Figure 4. In this figure we depict the effect of successively allow- 
ing extra parameters to vary. The dashed line shows the best fit 
model where only f!jj is allowed to vary, the dashed-dotted line 
where (Ami^a) vary, and the solid line where all ((Imi^Ai™) are 
free to vary. 

4.3.1. Allowing w ^ —1 

One further cosmological inference on which we inves- 
tigated the possible implications of residing in an under- 
dense part of the universe, was the equation of state of 
dark energy. This is characterized by the parameter w, 
which relates the pressure and density of dark energy via 
p = wpc 2 , such that p oc a~ 3 ( 1+w \ 

For dark energy that behaves as a cosmological con- 
stant, the value of w is -1. This val ue is supported by 
super nova observations, for example (jWood-Vasev et ahl 
12001 use the ESSENCE supernova s urvey to constrain 



the value of w to 1.0 7 ± 0.09 (see als olAstier et al 



Kowalski et al.ll2008D. (Although see iKessler et al 



2006 



2009; 



Sollerman et al.l I2009L for recent developments showing 
a potential deviation from w = —1.) Here we examine 
whether a void could significantly change the inferred 
value of w from the actual value. 

The LTB model we used as a description of our uni- 
verse was the same as that used in section 14.31 with 



M 



0.21 



QOUt 



0.3, n A = 0.7, and r = 70/i" 1 Mpc. 



In this LTB model the equation of state of dark energy 
is w = — 1, so if upon fitting models with varying w, we 
find a best fit different to w — — 1, we can say that this 
is solely due to the local under-density. 

At first, we simply add an extra parameter w, to our 
fitted models, so that each of Hm, ^a, and w are allowed 
to vary freely and independently. The error in w is very 
large, 44%, and the models fitted correspond very poorly 
to the external universe, the best fit parameters being 
(r2M,^A,w) = (0.314,0.53,-1.44). Such parameters are 
preferred since the low value of w leads to a sharp drop 
in distance modulus as z goes to zero. This mimics to 
some degree the effect of the void, depicted in figure |U 

Since fitting for variable w using only supernovae with 
no other constraints has large degeneracies, we now go 
on to demonstrate the behavior of the fitted value of w 
under certain prior constraints. 

First we adapt our fitting process to more closely corre- 
spond to conventional methodology. When determining 
the value of w, it is common practice to assume a flat uni- 
verse. We introduce the constraint Om + = 1, such 
that f^A is no longer allowed to vary freely. With this 
restriction the error in w is much reduced at only 7%. 
The constraint on fijyi and £l\ means that low values of 
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w can no longer combine with the required (^Mt^a) to 
mimic the void at low redshifts, hence this undesirable 
behavior is avoided and the errors in w are far smaller. 

Thus we find that using a prior from other observa- 
tions that constrains the wCDM model to be flat makes 
the result more robust to low-z density fluctuations than 
fitting a general wCDM model with no prior. However, 
one must also ensure that the prior used would not be 
rendered invalid by a putative local void. 

We note that the effect of a local void for the un- 
constrained and flat wCDM models is to push the best 
fit value of w down, and so could give the illusion of 
phantom-like dark energy (w < — 1) if the true value of 
w = — 1. 

Secondly we include a prior constraint on the value of 
Ojvl = 0.27±0.03. This is included in our fits by mod- 
ifying the x 2 estimate to xLtai = xIn(^m, ^a, w) + 
xUot&m), where ^ = (^) 2 . The best fit 
parameters under this constraint are (Qm,Qa,w) = 
(0.27,0.86,-0.85). Unlike the two previous varying-u; 
fitted models, the inferred w increases and no longer sug- 
gests a phantom like dark energy. 

5. THE LOW-REDSHIFT CUTOFF 

When fitting cosmological models to supernova data, 
it is common practice to remove data points of super- 
novae below a certain redshift. One reason is that at 
low redshifts the peculiar velocities of galaxies are a sig- 
nificant fraction of their recession velocities. Peculiar 
velocities therefore add a large amount of scatter about 
the magnitude- redshift relation at low redshifts. Typi- 
cally supernova studies have ne glected data from sources 
with a redshift below z n ~ 0.02 rtWood-Vasev et al.ll2007l : 
lAstier et alJl2006t iKessler et al.ll2009l ). since at this red- 
shift the mean peculiar velocities contribute less than 5% 
scatter about the mean recession velocity. 

A potentially more significant reason from the perspec- 
tive of the discussion in Section. [1] is that light from 
nearby sources does not give a good indication of the av- 
erage density of the homogeneous universe, since much 
of its path has been through the local Hubble-bubble, 
should one exist. This is potentially more significant 
than random peculiar velocities because it introduces a 
systematic shift onto all redshifts, not just a scatter. The 
value of redshift below which we neglect data is referred 
to as the low-z cutoff and we de note it by z . Figure 17 
of lJha. Riess. fc Kirshnerl ([200 7) shows the derived value 
for the equation of state of dark energy, w, can change by 
as much as 20% when the low-redshift cutoff is changed 
from zq ~ 0.008 to zq ~ 0.025 in the presence of a void 
that extends to the larger of those two redshifts. 

A commonly held notion is that by neglecting data 
from sources that could be situated within a Hubble- 
bubble, we remove the im pact of a local void (e.g. 
Una. Riess. fc Kirshnerl 120071 footnote 20). By neglect- 
ing nearby data the remaining light will have originated 
in, and primarily traveled through, the external region 
and thus should give a good indication of the external 
density. Here we investigate whether this low-z cutoff is 
sufficient to satisfactorily mitigate the distorting effect 
of the local Hubble-bubble. We want to know whether 
neglecting enough low-z data allows us to recover the cor- 
rect matter-density for the homogeneous region outside 



the void in our simulated data. 

We put this notion to the test by introducing a low-z 
cutoff, and gradually increasing it to see if the derived 
cosmological parameters converge to the input parame- 
ters outside the void. 

We note from the outset that even the light originat- 
ing from far away has to traverse the under-dense region 
before being observed, so removing nearby sources will 
not completely remove the effect of a void. 

Figure [5] shows four examples of different models that 
best fit the data as we move to progressively higher low-z 
cutoffs, z = [0.06,0.105,0.195,0.4] (Large-fi M case). It 
is clear that, as expected, the parameters of the best fit 
ACDM models are in better agreement with the asymp- 
totic values of the LTB model as we progressively neglect 
more low-z data. What may not be expected is that the 
best fit cosmological parameters do not converge to the 
known external density even after all data from within 
the Hubble-bubble have been rejected. 

We summarize this in Figure [6] which shows how the 
best fit f^M changes with zq. We see that as we remove 
data from sources within the Hubble-bubble, we do in- 
deed get a better indication of the external matter den- 
sity (which from our model we know to be fijjj* = 1.0). 
However, this figure also shows that the best fit J7m 
asymptotes to a value less than SIJJ 1 '. Regardless of how 
much data we remove we cannot deduce the true external 
density in this manner. 

We rep eat this analysis for the Small-A case. Fig- 



7(a)| shows that for this smaller under-density, the 



best fit mass density and cosmological constant density 
asymptote to the actual external densities as the low-z 
cutoff is increased. However, at the low redshift cutoff 
currently in use, z = 0.02, the error in Qm is 0.06%, 
and the error in J7a is 0.99%, the latter of which remains 
significant. 

An error of 0.06% in 17m is not likely to be a significant 
fraction of the observational error in the foreseeable fu- 
ture, but the error of 0.99% in 57a is already comparable 
to the observational error. This suggests that a higher 
low-z cutoff may be advisable. Figure [8] shows the low-z 
cutoff needed to reduce the errors in S1m and Qa to 0.1% 
of the exterior value. Assuming a void size of 70/i -1 Mpc 
(^lOOMpc), a low-z cutoff of about 0.035 Mpc is needed 
to reduce the error in J7a to 0.1%. Of course this is very 
dependent on the size of void we chose to test. We show 
in Fig. |9] how the best fit J7m and Q\ values change with 
the low-z cutoff for a range of void scales. 

Finally, Figure [TU] shows the best fit value of w as a 
function of zq. As was the case for J^m and the best 
fit value is different to the true value, but converges to 
the true value as we increase zq. When all three of 51m, 
Oa, and w are allowed to vary freely the error in w can 
be reduced from 44%, with no low-z cutoff, to 1.2% with 
a low-z cutoff of 0.035. At the currently preferred low-z 
cutoff of 0.02 there remains an error of 8.5% (again this 
is dependent on the void size chosen). When we impose 
flatness by setting 17m + ^a = 1-0 we effectively reduce 
the number of free parameters by one and as shown pre- 
viously, the error even with zq = is much reduced at 
6.7%. When we also add a low-z cutoff of 0.02 this is re- 
duced to a 1.4% error and a low-z cutoff of 0.035 reduces 
the error to just 0.2%. 
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Figure 5. These plots of distance modulus (normalized to the empty universe) show how the best fit model changes as the low redshift 
cutoff changes. Four models are plotted. In each case they are plotted as a solid line over the range of redshifts used in the fit, and then 
extrapolated as dotted lines over the range of redshifts that were excluded from the fits. As we implement progressively stronger low-z 
cutoffs it is clear how the best fit model deviates ever more strongly from the data in the nearby under-dense region. 



6. DISCUSSION 
6.1. Likelihood of Drawing Incorrect Conclusions 

We have seen throughout that the each time we allow 
a new parameter to vary, a better fit is found, but that 
the best fit parameters are altered such that they do not 
give an accurate representation of the cosmological pa- 
rameters outside the under-dense region. For example, 
in the matter only LTB model, allowing £7a to vary im- 
proves the fit but results in an f^A^O, despite the LTB 
model having SI a = 0. 

It is important to consider whether a cosmologist mak- 
ing such fits would consider an improvement in agree- 
ment with the data upon varying a new parameter, to 
be strong evidence for a value of that parameter differ- 
ent to that when it is not allowed to vary. Fitting an 
additional parameter, as opposed to setting it to a fixed 
value, will always improve the goodness of fit of a test 
model, regardless of whether this parameter is physically 
motivated. For example, a higher order polynomial will 
always give an equal or better fit than a linear function 
to any data set, even when the extra parameters are spu- 
rious. 

Clearly it does not give any extra insight to include 
extra parameters without any physical meaning. So to 
determine whether a parameter, or a particular value of 
a parameter does indeed well reflect the system under 



consideration, one can look at the improvement in the 
goodness of fit (GoF) resulting from varying the param- 
eter. Here we use the definition of GoF as x 2 /c?o/, where 
dof is the number of degrees of freedom (number of data 
points minus number of parameters). If the goodness of 
fit improves by a large amount upon allowing a param- 
eter to differ from a certain value, then this is taken as 
strong evidence that the underlying system is described 
by a parameter with value different to the previously set 
value. However, if the goodness of fit improves only a lit- 
tle, then it is likely that the improvement is simple due 
to the extra freedom, granted by the variation of the pa- 
rameter, for a model to fit the data set and its associated 
errors. 

This is often quantified in "Information Criteria" such 
as the Akaike information criterion, AIC = — 21n£ + 
2k, where C represents the likelihood, k represents the 
nu mber of param eters, and for gaussian errors — 21n£ = 
X 2 (|Akaikeill97l) . A A AIC > 6 between two models is 
considered significant evidence that any extra parameters 
in the model with the lower AIC are well justified by 
the improvement in the fit. In other words, it is unlikely 
that such a large improvement to the fit would occur 
unless the extra parameter described a genuine feature 
of the data. For a discussion of information criteria in 
a cosmological context and the more rigorous Bayesian 
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Figure 6. Upper: As we progressively remove data points from 
low-z, the estimation of the external matter density improves, but 
does not converge to the true value. Lower: ACDM models with 
cosmological constant give a good fit to LTB models, despite there 
being no cosmological constant in the LTB model. However, as we 
remove low-redshift data points, this spurious cosmological con- 
stant is reduced. 



evidence see lLiddld (pool : ILTddle et alJ pool ). 

To yield goodness of fit calculations using our simu- 
lated data which would match those made by a cosmol- 
ogist in an LTB universe, we fist add a normally dis- 
tributed random error with standard deviation 0.2 mag 
to the distance modulus at each redshift, this replicates 
the observational error a measurement would incur. As 
before we then assign an uncertainty to this measure- 
ment of distance modulus of 0.2. With this simulated 
data we can then calculate the best fit cosmological pa- 
rameters and the associated goodness of fit. To remove 
the dependence of the value of goodness of fit and best fit 
cosmological parameters on the particular set of random 
errors generated, we repeat this process 25 times with 
different sets of normally distributed random errors to 
get an average goodness of fit. 

We first consider the model in section I4.31 a small 
under-density embedded in a universe with both matter 
and cosmological constant. The best fit ACDM to this 
model is (0 M ,^a) = (0.305,0.72), with X 2 /dof=1.00. 
We note that the best fit matter and cosmological con- 
stan t den sities are slightly different to those stated in sec- 
tion I4.3I since here we have added many different sets of 
observational errors and taken an average, in section l4~3l 
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Figure 7. Same as Fig. [6] but for the realistic sized void consid- 
ered in Sect. I4.3I The effect is smaller, but even for this size void 
the error is not negligible given the expected precision of upcoming 
surveys. 
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Figure 8. This plot shows the low-z cutoff needed to reduce the 
error in and Qa to 0.1% for various sized voids. The more 
stringent constraint comes from Q\. 
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Figure 9. How the best fit values of Qm and Q\ change with low 
redshift cutoff for a range of void sizes. As the void size increases 
the low-redshift cutoff needed to reduce the systematic error be- 
low a certain threshold increases. The vertical lines represent the 
cutoffs needed to reduce the systematic error below 0.1% in each 
case. Void sizes tested range from 25 to 200 Mpc, and in each 
case the density contrast is taken to be 5 = —0.3. These are con- 
sidered amongst the plausible range of density fluctuations for our 
universe. 



we did not add observational errors, nor perform multiple 
runs. If we fit a set of models with only matter, 51m = 0, 
we find that the best fit CDM is (fi M , ^a) = (0.214, 0.0), 
with x 2 /dof=1.07. Here we see that allowing 51a to vary 
from 0, the goodness of fit is much improved (AAIC=19), 
and the cosmologist would conclude that there is strong 
evidence that the universe has a cosmological constant 
with SIat^O, which is indeed the case. This shows that 
a model with 51a = is unable to give a good fit to the 
LTB model with cosmological constant, and allowing 51a 
to vary is both required in order to obtain a good fit and 
does shed light on the nature of the actual universe i.e. 
predicts a non-zero cosmological constant. 

Now we give a counter example, where an improved 
fit should not be considered significant. In section 14.11 
we obtained the best fit ACDM model to a matter-only 
LTB model. With simulated errors, the average best fit 
parameters are (51m, 51a) = (0.85,0.41), where we have 
allowed 51a to vary from 0, despite the LTB model hav- 
ing S1 A = 0. For this ACDM , the xV^o/is 0.994 (AIC= 
301.2). Fitting a matter-only model (by setting 51a =0), 
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Figure 10. As with flm and f2j\> the under-density leads to an 
error in the best fit value of the parameter w, which is lessened as 
zq increases. When the model is constrained to be flat the error 
in the best fit w is smaller, primarily because the uncertainty in 
w becomes larger. Given the increased uncertainty in the curved 
case, the error is no more significant than the error when the model 
is constrained to be flat. 



we obtain best fit parameters of (Hm^a) = (0.87, 0.00), 
and a x 2 /dof of 1.009 (AIC= 304.7). As before, there is 
a clear improvement in the goodness of fit by allowing 
SI A to vary from 0, despite the fact that the model which 
generated the data had 51a = 0. This is a case in which, 
solely due to the presence of the local under-density, the 
astronomer may incorrectly conclude that there is evi- 
dence for a cosmological constant with density SIa^O. 
However, the AAIC between the two models is only 3.5, 
which should only be considered weak evidence in sup- 
port of the better fit. In this case the AIC indicates that 
there is a non-negligible chance that the extra parameter 
is spurious. 

Finally, we consider whether the presence of a local 
under-density would induce an astronomer to incorrectly 
conclude that the equation of state of dark energy was 
different to w = —1. Using the SmallA model as the 
data set, in section 14.31 we fitted a model where w is 
fixed as —1. Doing so again including observational 
errors we obtain a best fit universe of (51m, 51a, w) = 
(0.305,0.72,-1) and a { X 2 , X 2 / dof) of (299.0, 1.000). 
Allowing w to vary away from -1, the best fit model, 
as described in section 14.3.11 is now one with parameters 
(Q M , 51 A , w) = (0.246, 0.75, -1.65) and the \ 2 is reduced 
to 298.2, whilst the \ 2 1 dof basically stays constant (it 
worsens very slightly to 1.001). Despite the drastic de- 
viation from the w = — 1 model, the best fit is only a 
small improvement and AAIC= —3.2. Thus we would 
conclude that there is no evidence for an equation of 
state with w ^ —1, but rather evidence that the extra 
parameter gave too much freedom to the model leaving 
it poorly constrained. 

Thus we conclude that a void gives a high-risk of fool- 
ing us into believing an incorrect 51a but gives a low- 
risk of fooling us into believing io^-l. However, once 
we constrain to flatness we note that Flat-wCDM and 
ACDM have the same number of parameters and are 
equally good fits to the data, so we have no statistical 
reason to prefer one over the other. 

Table [2] summarizes the results presented in this sec- 
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Figure 11. When there is an under-density in the local universe 
all incoming data is rcdshifted. This is why removing the low-z 
data that come from inside the void does not allow us to completely 
recover the correct model. By effectively shifting all the points in 
the magnitude redshift diagram to higher redshift the void has 
changed the slope of the relation, which is the feature we're fitting. 
Therefore as we progressively remove data points from low z, the 
estimation of the external matter density improves, but does not 
converge to the true value. 

tion. They show the improvement in goodness of fit for 
each extra parameter added, which as discussed gives an 
indication of our susceptibility to be deceived by the void. 

6.2. Intuitive understanding of the void effect. 

The under-dense region has two effects on light pass- 
ing through it. Firstly, the lower density gives rise to a 
greater rate of expansion, H in > i? ut> because of less 
gravitating mass. Thus light traveling through this re- 
gion is cosmologically redshifted more than it would have 
been had it traveled through a medium of 51m = fijf'- ^ n 
addition, the under-dense region is at a higher gravita- 
tional potential than the external medium, so it causes an 
additional gravitational redshift to the light immediately 
before we observe it. Both these effects cause incoming 
light to be more redshifted than had it simply traveled 
through a homogeneous medium of density 51m = n^'. 
Hence the ACDM model we fit to this data will be one 
that redshifts light more than one with 51m = ^Jf' ■ This 
behavior is characteristic of a model with 51m < fi™' an d 
11a > 0. Thus by not accounting for a Hubble-bubble we 
mistakenly deduce a lower-than-actual matter density to- 
gether with the existence of a cosmological constant. 

In Fig.[TT]we show diagrammatically how a low-z cutoff 
is inherently incapable of removing the effect entirely. 
The extra redshift due to the void will be the same for 
all light sources situated outside the void. The effect is to 
shift the magnitude-redshift relation along in the positive 
z direction. In this figure we see that at each value of z 
the magnitude-redshift relation distorted in this way has 
a gradient more similar to that of a ACDM model with 
a lower CIm- 

This graphical interpretation of the effect is equivalent 
to realizing that even high redshift sources are affected by 
the under-dense region. This is apparent from the out- 
set, since light from far away still has to pass through the 
under-dense region in order to reach us and be observed, 
and hence will incur a greater than expected redshift, 
if we do not take the local inhomogeneity into account. 
Thus, it is evident that the current practice of neglect- 



ing sources at low redshifts in the hope of removing the 
effect of the void is optimistic, but as we have shown 
quantitatively, this approach does have the desired effect 
provided that the under-density is of the scale of ob- 
served inhomogeneity in our universe. For larger voids, 
removing supernova observations from low redshifts will 
not remove the effect of the void, regardless of how much 
data we neglect, as demonstrated in sections |4~T1 and l4~2l 

Although we have only simulated under-densities in 
this paper the same principle applies to over-densities. 

We have also tested whether the asymptotic value of 
51m relates to any average matter density in the volume 
we can define, including the lower density of the void. 
The answer is no. 

7. CONCLUSION 

We have shown that a local density fluctuation can 
leave a significant systematic error on our cosmological 
inferences from distance measures. This applies not only 
to supernovae, as discussed here, but also to measures of 
angular diameter distance such as baryon acoustic oscil- 
lations and the cosmic microwave background. The low-z 
cutoff imposed on supernova data is absolutely crucial to 
avoid significant systematic errors, far above the statis- 
tical uncertainty in current observations. This is par- 
ticularly true for the cosmological constant, whose form 
can mimic that of a local void. We have shown that 
for an under-density with a scale and depth typical of 
voids in the current standard model and the low-z cut- 
off currently in use, z=0.02, a systematic error of 0.99% 
remains in the deduced value of 51a and hence a slightly 
higher low-z cutoff of 0.035 is advisable, at which the er- 
ror drops to a negligible 0.1%. The equation of state of 
dark energy is less susceptible to bias induced by a void 
than a cosmological constant, because it does not easily 
mimic the void effect, but a local under-density can still 
induce wild discrepancies in the best-fit values. 

We have also shown that introducing a low-z cutoff 
can reduce but not remove the systematic error on the 
inferred cosmological parameters. Thus a void imprints 
an irreducible error that can only be removed by fitting 
the extra parameters needed to allow for a local density 
fluctuation. The magnitude of this bias is small com- 
pared to the current observational uncertainty on 51m, 
51a, and w. Nevertheless, as our cosmological inferences 
become more precise, this bias will no longer be negligible 
and will need to be accounted for. 
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